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There are investigated the problems of the flow of a supersonic 
jet out of a vessel with plane side walls and the problem of the super- 
sonic flow about a wedge when there is a zone of local subsonic veloc- 
ities ahead of the wedge. 


INTRODUCTION 

In the present paper it is assumed that the reader is acquainted 
with the work of S. A. Chaplygin ("On Gas Jets" (reference 1)) and with 
the method of computation of plane -parallel supersonic flows given by 
Prandtl and Busemann (reference 2, see also references 5 and 6). There 
ia recommended a preliminary acquaintance with the work of F. Tricomi 
"On second order partial differential equations of mixed type" (refer- 
ence 3) whose methods undoubtedly will be capable of being used in 
proving the existence of the solution of the problems considered by us. 

Since in what follows we shall everywhere make use of the notation 
of Chaplygin we shall here present the formulas and notation of impor- 
tance to use. Chaplygin makes use of the method of the holograph. As 
the independent variables he chooses in the first place the nagnitude 



where V is the flow velocity at a given point, V m is the maximum 
velocity corresponding to the stagnation temperature Tq (that is, 
the temperature arising in front of an obstacle in the flow) char- 
acteristic for the given flow; V m is given by 

V 2 m = 2Jgc p T 0 (2) 

where J is Joule's constant, g the acceleration of gravity, Cp 
the specific heat for constant pressure, and Tq the absolute 
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stagnation, temperature. The second independent variable is taken to 
be the angle of inclination of the velocity 8, Written in these 
independent variables the stream function '4 / in the case of irrota- 
tiona l flow satisfies the equation - 

5 f 2 t gj/1 1 - (20 + 1) T a 2 ^ n /,x 

[(1 - t) 13 ^ T j 2 (1 - t) (3+1 S0 2 

where 

e ■ i?”r < 3a ) 

and 

K = “■ (3b) 

G v 

the ratio of specific heats at constant pressure and constant volume , 
respectively. 

The value 


- ~ U 1 


(4) 


corresponds to the critical velocity; that is, the velocity of the 
flow equal to the corresponding local sound velocity. 


On introducing the auxiliary variable 


C = 


n (2p+l) -1 . 


(1 ^t) 

2T 


dT 


equation (3) assumes the form 


3a 2 


3^ 

+ K —— = 0 

3e 2 


( 5 ) 


( 6 ) 


where 


K 


1 - (20+1) t 


(1 " T) 


23+1 


(7) 
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Thus, equation (6) for subsonic velocities will "be of the elliptic 
type and for supersonic velocities of the hyperbolic type. 

Chaplygin further considers particular solutions of equation (3) 


of the form 



where 

^y (t,' 0) a Z y (t) sin 2v6 

(8) 


z v (t) = T v yy (T) 

(9) 

and yy ( t ) 

is the hypergeometric function 



yy (t) = F (ay, by; 2v + 1; T ) 

(10) 


where 


+ by ~ 2u ~ 3, ayby — — 3"^ (2u + l) (10a) 

In the theory of Chaplygin an important part is also played by 
the auxiliary function Xy (t) 




1 + 


t y' 

v yp 


u 


i fin 

u Z„ 


( 11 ) 


The problems considered by Chaplygin for flow velocities remaining 
everywhere below the velocity of sound reduce to the problem of 
Dirichlet and are solved with the aid of series combined from the 
special solutions of the form (8). To what boundar y problems for 
e qua t .1 o n (3) the pr ob lem s of Cha plyg in reduc e for mi xafl— &uh~- .and. supe r - 
sonic flows r em ained u nknown] "Basing himself on the work of Tricomi 
(reference h), the author has succeeded in finding a formulation of 
these problems and to e stabli sh the uni queness of their solutions. 

In what follows the author hopes to give a mathematically well 
founded and practically suitable solution of the problems stated. 
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I. SEDUCTION OF THE PBOBLEMS OF THE FLOW OF A SUPEBSONIC 


.1ST TO THE PBGBEBM OF TEICOMI FOP THE EQUATION OF 


CHAPLYGIN UNIQUENESS THEOBEM FOE THESE PBOBLEMS 


The problem of Tricomi is bhe following: Let there by given a 

linear partial differential equation of the second order which on one 
side of the curve C in the plane of the independent variables is of 
the elliptic type and on the other side is of the hyperbolic type.. 

Let us consider the region D bounded by the curve L lying in the 
elliptic region with its ends lying on the curve C and with the 
characteristics x-j_ and Xg belonging to different families and 
starting from the ends of the curve L (fig. l), Let the values of 
the solution be given on the curves L and X] but not on Xg . 

There is sought a solution in the region D. 


This boundary problem was first formulated bjr F. Tricomi (refer- 
ence 3) as applied to the equation 



& 


o 


(i) 


Tricomi proved the uniqueness and existence of the solution of 
this problem. In this section we shall reduce the problem of the 
flow of a supersonic jet to a certain problem of Tricomi for the 
equation of Chaplygin (see introduction equation (6)): 



Zd 2 


( 2 ) 


The coefficient K for small T or large o is equal to unity 
and drops with increasing T (decreasing a). For t = (l + 2p) - ^- 
(or a = 0) we have K = 0, and for T > (1 + 2p) -1 (or O < Q) 

K< 0, Thus equation (2) is of the elliptic type for a > 0 and 
hyperbolic for c < 0. We shall first prove the uniqueness of the 
solution of the problem of Tricomi for equation (2). 

We consider in the plane (0, o) a finite region bounded by the 

curve ABC lying in the half plane a > 0 and the characteristics 
AD, CD, lying in the half plane o < 0 (fig. 2). We assume that 

in this region the solutio n \l' of equation (2) is taken equal to 

zero on ABC ancT'on "CTb We shall show that this solution is equal 
‘ro'zero over" the 1 entire ''regfion. 
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We consider first the solution in the triangle and will show 

that 



> 0 

0=0 


(3) 


We transform the equation 



hy integrating hy parts. We have 



> (5) 



Along the characteristic DA we have 

de - J | k| do = o 

> 


do _ d0 = J \ K I do 

J |K| 


( 7 ) 
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Hence 


/'*'(! ae + |k i H to ) - (I ao + 

DA DA 


io) = 
de / 


J' \/Jk~i d'l = 



= o 


o 

Cmin 


DA 


so that 




> 0=0 


d© 


1 

2 


t 

J 


DA 


5JH* 2 ao 

da 



We now compute : 

dv / i~K| _ 1 d |Kl /_ 2r \ = T 2p (2p + l) ■ 

da " ZjW\ dT \ (1 - T) P / (1 - tJMkF (1 - r) 2i3+2 


Thus, to prove the inequality (3) it remains to show that 



d© dO ^ 0 


To prove this we shall transform to characteristic coordinates 

dA = d© + a/[h] da*) 
dp. = d© — VTST da f 


( 8 ) 


d© da 

O) 


■ < 0 
( 10 ) 


( 11 ) 


(12) 


from which we obtain the Jacobian determinant: 
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Hfc-gj- * - 


Differential equation (3) "becomes 

n 2 


_ 4 | K| JOL + = o 

§A 5|a da \SA i/ 


On the other hand 




\ae, 


aA Bn 


or 


M - |K| 


W 6 

v C> 0 


do da = - 2 


JJ 


— -r-a^ . 

K| 3A5: dX ^ 


where the integral on the right is taken over the area CA'D 1 
To compute this integral we rewrite equation (14) in the 

form: 


where 


/TFT (^L _ 

v 1 1 w a u 


M (a) 


a 2 ^ 

a a an 


M 


(a) .41* 


,3/2 


d J | K ! 
da 


Then 


JJ ^ ll S5 dA dd -ff A* r (|J dX dd *JJ M I 


a^ a 2 ^ 


an aA a^ 


dA 



i 

J 


"m ** 



A=n 


dn 


A=0 




(13) 

(14) 

(15) 

(15a) 

(fig. 3) 
following 

(16) 
(16a) 

dn = 

( 17 ) 
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since for small o 


M = 0 (o' 5 ) 


and for continuous 6 , b'ty/cMJ 


a* 


o 0or l/2 ) 


(17a) 


(17b) 


On the other hand. 


</|K| - 


S - sm ^ 


(2 + p) T - 2 


4-v/fKl uu 3 (2P + 1) t 2 

This expression is negative if over the entire triangle ACD 

2 


(18) 


T < 


2 + p 


(19) 


The above inequality expresses the fact that the Mach number M should 

2 

be less than 2: for, with T = = 

2 + 6 


M 2 = 


K - 1 1 


2 2 (K - 1) 


K - 1 2K - 1 


~ (2k - 1) = 4 (19a) 


This means for K = 1.4 the base 9q of the triangle must satisfy 
the inequality 

8 0 <SS 54° (19b) 

Whether this restricting condition for the proof of the inequality (3) 
is essentially required or whether it is only connected with our 
method of proof is as yet unclarified*. 

Let us now consider the region ABC. By integrating by parts 
we obtain as above : 


*The proof remains valid for any equation of the form (2) where 
K is a regular function of o for a = 0 and dK/do > 0 for o = 0 
and K(0) = 0, provided that a in the triangle ADC remains suffi- 
ciently small. As applied to equation (l) the proof remains valid for 
any size triangles ADC. The same is true of the proof of the unique- 
ness as a whole. 
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and taking account of inequality (3) 


( 20 ) 


( 21 ) 




|a=o 


de 


o 


( 22 ) 


From equations ( 24 ), (22), (11), (17a), and (19) it follows now that 


\|/ = 0 


(23) 


as was to he proved. 

The proof of uniqueness here given is applicable under the con- 
dition that the transformations encountered are valid. This com- 
pletion of the uniqueness proof we shall give after investigating the 
properties of the solution of our problem. 


Let us return to the problem of the flow of a supersonic stream. 
We consider a vessel with symmetrically arranged walls forming an 
angle 20° (fig. 4(a)). We assert that for a sufficiently small 
external pressure the problem is reduced to the following problem of 
Tricomi: In the region OA'B'C (fig. 4) a solution ^ is sought 

under the conditions: 


1 ° 

2 ° 




\J/ = - ^ on OA'B' 


> 


\j/ = 0 on 0C 


( 24 ) 
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Similarly on tlie "boundary of the region OABC we must have 


1° ii - | on CAB j 

\ 

2° ^ = 0 on 0C J 


(24a) 


The obtained solution gives a mapping of the region OA'B'CBAO 
in the plane (x, y); thereby are obtained the region of subsonic 
velocities, the curve of sound velocities, and the part of the super- 
sonic stream touching this curve. The continuation of the supersonic 
stream to infinity is obtained by the method of Prandtl and Busemann 
(reference 2) . 


We now proceed to prove the above statement. According to the 
conditions (24) and (24a) the walls of the vessel correspond to the 
radii 0A' and 0A, and the axis of symmetry of the vessel to the 
radius 0C. As regards the characteristics AB and A'B’ there 
correspond to them in the (x, y) plane the points A and A‘, the 
opening edges of the vessel. For, if along the characteristic of 
equation (4) the stream function \j/ is constant, the pobencial c? 
is likewise constant (reference 4, section I, formulas 1.15 and 1.16). 
But if along a certain line cp and are constant, then the 
coordinates x and y are also constant. It remains to show that 
the obtained flow may be continued in the form of a stream with con- 
stant pressure (constant velocity) on its "boundaries. Such contin- 
uation is possible if the ratio of pressures Pi/pq is less than 
(not equal to) a certain function of the angle 6: 


Pi 

P 0 



(%) 


(Pg — pressure at point B) 


(25) 


We give below a table of values of this function: 

0 Q = 10° 20° 30° 40° 50° 54° 

f (9 0 ) = 0,33 0,26 0,20 0.17 0,14 0,13 

Let us consider first the case p^_ = p-g. We draw the arc of a 
circle BB' with center at the origin of coordinates, and prolong 
the characteristics A'B’ and AB to their intersection D, 

We now find the solution V = °1 "the equation of Chaplygin 

in the triangle B'CB passing along the characteristic CB into the 
previous solution ^ and equal to (+ Q/2) elorg BB' . Further, 
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ve find the solution \j/ = in the triangle BB'C passing into the 
solution ^ along CB' and equal to (- Q/2) along B'B. 

Further , ve find the solution ^ = ^4 in the quadrilateral CB’DB 
passing along CB inbo ^3 and along CB' into ^2- Continuing, 
ve find the solution ^5 equal to \|/ 4 along B'D and to + Q/2 
along B'B, and symmetrical to the latter, the solution ^g equal to 
\|/ 4 along BD and equal to (— Q/2) along BB' . We then find the 
solution ^7 equal to ^5 along DB and equal to tyg along DB' , 
etc. 


The regions in the plane (x, y) corresponding to these solutions 
are denoted in figure 4 by the corresponding numbers. Thus ve evidently 
obtain the flow with pressure pg on the boundary. 

If the pressure in the outer region is less than pg we proceed 
as follows (fig. 5). We draw the arc of a circle EE' with radius 
corresponding to the pressure Pp < Pg« The points E and E' must 
lie on the prolongations of the characteristics AB and A'B'. The 
intersections of this circle with the prolongations of the character- 
istics CB and CB' we denote by F and F' , We draw finally 
through F and F' the conjugates of the characteristics intersecting 
in the point D on the axis of u. 

We now find the solution ^2 i- n the quadrilateral CHEF' that 
passes into along CB and is equal to Q/2 along BDF ’ ; then 

the solution ^3 in the quadrilateral CB'EF that passes into 
along CB' and is equal to (-Q/ 2) along B'E'F'. Further, we find 
the solution \|/ 4 in the qualdri lateral CF'DF which passes into ^2 
along CF' and into ^5 along CF. There are then found two 
solutions in the triangle FF'C, etc. as was shown in the previous 
case. We thus obtain the flow with the constant pressure Pp < Pj$ 
on the boundary. 

To prove the existence of a steady continuous supersonic stream 
flowing out of a vessel it is necessary only to prove the existence of 
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the solution of the problem of Tricomi*. A strict proof of existence, 
as ha3 already been said in the introduction, has not yet been obtained 
by us. The fact, however, that the solution of the problem of Tricomi 
for equation (i) exists and for equation ( 2 ) there has been shown the 
uniqueness of the solution makes it appear probable that the solution 
of the problem of Tricomi for equation (2) likewise exists. 

It should be noted, however, that the proof of uniqueness was 
obtained only for the values 9 q $ 54°. If this corresponds to the 
actual state of affairs and if the existence of the solution were 
established only for the values 0 q 5 ? 54° this would mean that the 
flow out of a symmetrical infinite vessel with straight walls is 
possible in the form of a steady continuous supersonic flow provided 
these walls include an angle not larger than 108°. The assumption is 
here made that for 29 q > 108° supersonic flow without density jumps 
(shock waves) is impossible. It would be interesting to check this 
assumption experimentally. 

With regard to the obtained solutions the curves of the velocity 
of sound start from the edges of the opening. It is to be noted that 
for p = pp there corresponds to the characteristics A’B ] and AB 
one point of the plane (x, y), namely, the edge of the opening and 
this is also true for the case p-j < p^ with the corresponding char- 
acteristics A'E' and AE . This means that the flow in the neighbor- 
hood of the edges of the opening has the character of a Prandtl -Meyer 
flow (reference 5), that is, the character of the flow about a corner 
with expansion. The angle of inclination of the boundary of the jet 
as compared with the direction of the wall should be not less than 

%/ 2 - 


The flow within the vessel, since it is entirely determined by 
the solution \|/p of the problem of Tricomi, does not depend on the 
outside pressure pp provided Pp Pg« Hence the quantity of air 


This proof must of course be completed with the proof that the 
Jacobian D(x, y)/ D(u, v) or the magnitude (d^/cio ) 2 + K(5^/de ) 2 
for each of the solutions ^p, . . , has a constant sign. Other- 

wise the components would not be unique functions of the coordinates. 
In this case there would be expected the appearance of density jumps 
in the flow. It is not difficult to show to which types, according to 
Christianovich (reference 4) the flows considered belong. The flow 
\|/p in its supersonic part and also the flows and ^3 are mixed 

flows, the flow \[/p is a flow of rarefaction, the flows vj/g and tyg 
are mixed flows, the flow \J / 7 is a flow of compression, etc. 
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per second likewise does not depend on p-j_ as entirely corresponds 
to the well known experimental facts. 

By the above indicated methods it will not, however, be possible 
to find a solution if 



(26) 


In this case the problem is reduced to a boundary problem which 
is a generalization of the problem of Tricomi. The solution is 
sought in the region OCD'B'A 1 , where A'B’ and CD* are arcs of 
the ellipsoid of Busemann and B'D* an are of a circle corresponding 
to the given external pressure. The points 0, A', C* are the same 
as in figure 4(a). The boundary conditions are the following: 


\j/ = - | on OA'D'B' 


\J/ = 0 on OC 


(27) 


The uniqueness of our solution has not yet been proven but is 

K 

very probable. In the limit for p = p* = Po P* 1 ® a1aove 

boundary problem goes over into the Dirichlet problem and its solution 
into the solution of Chaplygin. 


II. REDUCTION <SF THE PROBLEM OF A SUPERSONIC FLOW ABOUT A WEDGE 
IN THE CASE CF THE FORMATION OF A SUBSONIC ZONE AHEAD OF THE 
WEDGE . TO A BOUNDARY PROBLEM FOR THE EQUATION OF CHAPLYGIN 
IN AN INITIALLY KNOWN REGION OF THE VELOCITY PLANE 

THE' OREM OF UNIQUENESS FOR THIS PROBLEM - ' 

P " 

In the case ' here considered the entropy behind the density, , jump 
is variable. In connection with this in the equation of Chaplygin 
for the flow -there appears a part on the right side proportional to 
the derivative of the entropy with respect to the stream function 


14 


MCA TM No. 1155 


(reference 5). In what follows we shall neglect this right part, as 
also in general the variability of the entropy T" The flow then remains 
potential (reference l) and the equations of Chaplygin remain in force 


3_cp 

30 3<j 

3cp = ^ 3^ 
3a - 30 


(l) 


3 2 ^ „ 3 2 ^ 

o + K r— 

3o 2 


0 


( 2 ) 


The problem of the flow about the wedge can now be reduced to a 
boundary problem in the region OAEDE of the plane (u, v) (fig. 6 ). 

In this figure OA is a segment of the u axis, AB is the arc of 
the strophoid giving the -velocity behind the wave front lying within 
the circle of the subsonic velocities. The equation of this strophoid 
(reference 5) is 




2 

K + 1 


Vp + 



u 


(3) 


Bf) and ED are arcs of the characteristics (epicycloid), B 
and E lying on the circle of sound velocity, and OE is the radius 
making angle 0 q where the latter i 3 the angle between the sides of 
the wedge in the direction of the approaching flow. The boundary con- 
ditions are : 


f = 0 on AOED 

(4) 

u f = V-g at pome B 

(5) 


where ^3 is assumed given. 

On the arc of the strophoid there must be satisfied such con- 
dition as would assure a continuous change in the stream function on 
passing through the wave front. 

The corresponding transformation of the (u, v) plane into the 
(x, y) plane is shown in figure 6 . As in the previous problem, to the 
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characteristic ID corresponds a single point E-D in the (x, y) 
plane. At this point _( corner at the hase of the wedge) there arise 
flows of the Prandtl -Busemann bype. The continuation of the flow 
beyond the Mach line OB is not of interest since it has no effect 
on the flow in front of this line (this continuation may be deter- 
mined by the method of Prandtl -Busemann) . 

/ 

Under these conditions there will evidently be satisfied in 
the (u, v) plane those boundary conditions in the (x, y) plane, 
which are a consequence of the formulation of the problem of 
Chaplygin*. The value of is proportional to the height of 

the wedge. It remains to render more precise the boundary con- 
ditions on the arc of the strophoid AB. 

Let p-[_ be the density in the undisturbed flow and X the 
angle of inclination of the density .jump at an arbitrary point 
(fig. 7). Let Vp be the velocity of the undisturbed flow, pQ 
the density at a stagnation point. We recall that pQd^ gives the 
difference in the discharge at two infinitely near points. Then 
along the discontinuity 

P 0 d^ = PpVpdy (6) 

On the other hand (reference 1) 

*-& i * + & w -T 29 ?- a * + 


sin 6 f cKp . , \ 

— 1,55 30 + Se a 7 ’ 


Po cos Q j,i, sin 9 ( Tr 

== — dv + — == — IK r— do — c— dS 

p v v \ d© do 


From equations (6) and (7) we have 


„ Sine U » do _ a as 

) p Q v \ be do , 


Further (fig. 8) 


(7) 


(8) 


-Jf 

Of the conditions which are satisfied on the wave front we 
have rejected one. This, however, is unavoidable since we have 
neglected the change of entropy. 
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Y 1 cos A = V cos (A - d) = Y s (9) 

(that is, the tangential velocities do not change on passing through 
the discontinuity) and 

PlVl sin A = pY sin (A — 0) (10) 


(that is, the flow discharge does not change in passing through the 
discontinuity) (reference l). Hence 


and 


1 cos 0 1 / 

PlVl pv" = Pl v x V 1 


s in ( A — 9) cos 0 
sin A 




Y 

sin 0 



cos 0 


( 11 ) 


JL cos A ^ __ sin (A — 9) cos 6 

sin 0 cos~(A ~~07 sin A 


ctg A 
Pt 


( 12 ) 


Hence the boundary condition (8) becomes 


< t- 


n 


<$Ai Bv PQ i 

K da — — - d0 = — ctg A d\|/ 
c >0 80 P x 


(13) 


Since along the strophoid 0 is a known function of a the equa- 
tion (13) gives a homogeneous linear relation between B^/B0 and 
B'l'/Bcn 


jj We shall now show that the conditions (4), (5), and (13) deter- 

/mine the solution of equation (2) in the re gion O AHDE uniquely , or 

in other words that the homogeneous conditions (4) and (13) determine 
the stream function except for a constant factor. For this purpose 
it is necessary and sufficient to show tha + the solution of equa- 
tion (,2) satisfying conditions (4), (13), and (5), '4(E) = 0 must 

be identically equal to zero. To prove thxS it is sufficient to show 
that from the satisfying of the condition (4) along AOE, (13) along 
AB, and (5) at the point B there follows 
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where the integral is taken along the line c = 0. For, in section 1 
it has already been shown that due to' the satisfying of the condi- 
tion \1/ = 0 along ED* 


'E 


'B 


\l/ 21 d© 5? o 
* Sa 


(15) 


and that from (14) and (15) we have 

\|/ s 0 

We shall prove inequality (14). We have: 


(16) 


o = 

J J \5o 2 bey 

QAEB 


r 


&0 do 


© 


' . K 


W 


d© da + 


- f * [ K I? 30 - 1 de ) * 


s^y 

day 


+ K 5 ^ 

\O0/ 


d9 da + 


r 'E 


5'4 


Ctg A • \|/ d^ - '4/ d0 


5a 


(17) 


AB 


'B 


Hence 


£* g “--JJ I 


+ K V5© 


d© da + 


p 0 

2p l 


~B 


J A 


\j/ 2 — — «c o 

sin 2 h 


(18) 


It is here necessary of course, as in section 1, to assume that 
at the point D M < 2 (or that the central angle of the arc B is 
less than 54°. 
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vhich proves the uniqueness theorem*. 

Since there was here assumed the validity of the transformations 
given it is necessary to supplement this proof of uniqueness by a 
proof of existence. To the problem of existence and the effective 
method of finding of a solution we hope to return later. 

If the pressure behind the wedge is greater than the pressure 
at the point D (fig. 6) the region in the plane of the holograph 
and the boundary conditions vary in the same maimer as in the case 
of the flow out of a vessel (see remark to section I). 


III. TWO LEMMAS TO THE THEORY OF THE EQUATION OF CHAPLYGIN 

POSSIBILITY OF APPLICATION OF SERIES OF THE 

TYPE OF CHAPLYGIN TO THE PROBLEM OF . 

A SUPERSONIC FLOW FROM A VESSEL 

In this section we shall prove two lemmas to the equation of 
Chaplygin which we intend to use later for a proof of the existence 
of the solution of the problem of Tricorn! . The first of these 
lemmas is an asymptotic formula for the logarithmic derivative of 
the function z v (t) of Chaplygin (see introduction, formula (ll)) 

for T = (2p + l)" 1 and large X> , namely 

z 'v (jjTT T 
(pTTT 

where C is a constant independent of u and the symbol 0 (l) 
means a bounded magnitude . 


1 = ct 2/3 


0 (1) 


( 1 ) 


*The proof of uniqueness is applicable only if the central angle 
of the arc BE is less than 54 . In this connection there is also 
obtained the restriction for the angle S L) (fig. 6), The limiting 
angle depends on the Mach number of the approaching stream. 

For M = 1, 9 0max = 54°, for M = oo ; ©Omax = ^ s °- question 

as to whether this limiting angle has physical significance still 
remains open. 
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In preparing the paper "The Theory ~of ~ the Laval Nozzle" for 
publication we were able to render this formula more accurate. We 
obtained an asymptotic formula for large b : 


'b 


= c // 3 + c 0 + Cl b - 2 / 3 + 


T = 


2p+l 


_2k 

+ C k b 6 + 0 



ov. o \ 



In the case t = 


m +“ i) 

inequality proven by Chaplygin, 


this formula is involved in an 
namely, the inequality 


/(Jl - 2(3s 


+ 23 s 


/f/p (1 + 2P) 2 

V 2b 2 


> > 



. ps 2 (1 + 23 ) 
+ _____ 


( 2 ) 


where 


(2a) 
(2b) 

We now proceed to the proof of (l). By definition of the func- 
tion z v (t), the function 

= z v (t) Bin 2b 0 (3) 


*b = 


1 + 1^_ 


z ’ 

i Li. 

b z-,, 


and 


s = l 


satisfies the equation of Chaplygin. If we replace the variable T 
by the variable 0 (reference 1, section V, formula (Si)) then from 
the equation 


to 2 


f K ~~ = 0 

de 2 


(4) 


there follows 

r v (a) - 4 b 2 K^(a) = o ( 5 ) 


where 


t’j(a) = z v (t) 


(5a) 
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We recall that 

Zu(t) = T v y u (T) (6) 

where y-^T) is a solution of the hypergeometric equation 

T(l - t) j" v + [2v + 1 + (p - 2y - 1) t] y' v + 0l?(2u + 1) j v = 0 

(7) 


regular for T = 0. 


Equation (7) has a second independent solution of the form 


y (2) v ( t) = t- 2 %(t) (8) 

where g^T) is regular for t = 0. Therefore equation (5) has a 
second independent solution 

i (z, v (a) - 2 (2 ) „(t) - T' v g„(T) (9) 


From the formula for the coefficient K (section 1, formula (4a)) 
it follows that near 0=0 

K = ao + ho 2 + . . . (10) 


where 


a = 2 



(10a) 


Since for T < (2p + l) ^ (or a > 0) K is hounded, it follows from 
equation (10) that 


Ik 


- ao | < Bo 


2 


for a > 0 


( 11 ) 


In differential equation (5) we now replace the coefficient K by 
its approximate value equal to aO. We obtain the equation 

l" v - 4b 2 aoL v = 0 (12) 

where £y(c) is that solution of equation (12) which for a-* 03 (T— »0) 
approaches zero and which is equal to t v for 0=0. This solution 
has the following form (reference 3, section III, formula (12)): 
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l v M 


\Wo) 


ACoT 


where 


l £r> 1_3 
-^ep-gp 


Ml) - I . !,r ir oo. (f + ^ |p) dp (K 

-'o \ / 

For what follows it is of importance that the function A(t) 
for any positive m satisfy the inequality 

N£)l < nr « > o) (li 

¥e now denote ty 5^ the function 

= ^u " ty (l* 

This function satisfies the nonhcmogeneous differential equation 

(sL)" -4v 2 E5tp = 4u 2 f(c) (!' 


where 


f (a) = (ao - K) £„(a) 


(17a) 


According to the general theory of homogeneous linear differ- 
encial equations there follove from (17) 

dL v (a) = if |t v (o) J° f(c') t (2) y (a') da’ + 

+ t^ 2 \,(a) r f(o') t v (o') do'V + c 1 tp(c) + g 2 ^ 2 \{o) (18) 
Jo J 


where 


r 

± v ^ v 

t W 

i»u b u 
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It is easy to prove that 


C 2 = o 

In fact, for small T (large a) 

£ v (a) = o (t u ), = o (t" u ), da = - — o (i) 

f (a) = o (a 2 ) = o ( in 2 t ) = o (t _£ ) 
where £ is an arbitrarily small quantity. Thus 


J 


f(o’) S (2) »(a') da' = 


, r 


(1+BP)-1 


-J. 


0 (t'-”- 1 -*) dr- 


J 


a 


f(o') 4(a’) d o' = [ o (x^ 1 ) da' = o (x u_€ ) 

J o 

From equations (16), (21), and (22) we have 

sL(c) = 0 (T' e ) + Q 2 l {2 \(o) 


Hence if C 2 4 0 then 

lim 5^p(a) = CO 

(J->c o 

which contradicts the definition of this function 
We shall now compute Cj_. We have 

f~ CO 

b£ v (o) = ~ ^ 2 \(o) I f(a') t v (a') da' + c 1 S v (o) = 

u o 


whence we obtain for the constant: 

f( 2) 

Ci = - 


4V“ 

A 




(o) 


f (o' ) L,(a') da’ 


o 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(23a) 

0 (24) 

(25) 
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and for the derivative 

st'v(o) = '£§-< t ( 2 ) 'y(o) J f (o') £ v (o') da' j> + Cjt^o) = 


Yjoj f(o') Sv(a') da' 


1T(oT 


nco 


Jo 




(a a' - e) £ v (a* ) da’ 


(26) 


so that 


t\(0) 

17 s i ’ 


$ 


4aV 


A’(0) 

TCcT 


4ir 

aToT 



Uo') 

'TJoJ 


da 1 


(27) 


However, according to Chaplygin (reference 1, vol. II, p. 30) 


0 < 


Uo) 


Uo) 


- ^ 1 for 0 < T ^ (1 + 2p) 


-1 


(28) 


(or 0 < a) 


lienee 




A(t) 


B£ 2 d £ 

(4a U 2) 2 7 5 ‘ (4aV 2)l/3 



( 29 ) 


Thus according to (27) 
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which proves formula (l). 

The computation of the magnitude X'(0)/A(0) from formula (14) 
and the known properties of F functions gives 


X'(0) 

TToJ ~ 




(30a) 


The second lemma is a consequence of the first and may he stated 
as follows: Let ^ be a hounded solution of equation (4) defined in 
the region 


O > C, 0 < 0 iC Q c (31) 

Let the limiting values for 3 = 0 and 6 = 0q he 

\!/(a, o) = \J/(a, Q 0 ) = 0 (32) 

Then there exists a kernel K( 0 i O') net depending on V the prop- 
erties of which are determined by the equation 

K( 9, Q') = A(l 9 - 9' | - 1 / 3 - (0 + O' ) _1 / 3 - (0 + 0' - 20 o )‘ 1 / 3 ) + 0(l) 


(33) 

which permits expressing the boundary values of ^ on the arc a = 0 
in terms of the boundary values of d^/do for 0=0 since 

r e o 

^(0, 0) = j K(e, O') vo, 01 ) d0, (34) 


To satisfy equation (34) it is required merely that the square of the 
function ^(O, 0) be integrable. 


We shall now prove this. We compute first the kernel K(0, 0’ ) 
assuming that it exists. In the particular case 


"it - % = T-JoJ 3in 2ve 


(35) 


we have 
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C^ v 

bo 


0=0 


M°> 


sin 2D0 


¥e introduce the notation 


v ~ 29, 


- e 
o 


X„ = 


r„(o) 


n ' 5„(o) 

Equation (36) may then he written in the form 


\k VT (0. 9) 0 

V 0 ' e) = — % = K ( 0 ; ^’)% a (0, 0’) as* 

n Jo 


or 


29 r 


K(0. O') sin 2n-i5'di 


•-'0 


sin 2n£' 

X n ' 


whence 


K(0 ; 0') = 


2 \ sin 2 nfl sin 2nd* 

0 O Z_, 


n-1 


'n 


It remains to investigate the convergence of this series 
According to equation (30) 


X n - 


1 

*n 


X'(0) 

X(o) 




4aU L 


x(0) 

X 1 ( 0 ) V' 4au 2 


[l + 0(n~ 2/3 )J 
[l + 0(n" 2 - /3 )J 


(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 
(42a) 


whence 
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K(e, o') 


Mo) 


CO 

V 


o A’(0)3/4a ^ 


> - 


oin 2nd sin 2nd' 


v 


ITT 


+ / sm 

n^T 


n=l 

2nd sin 2nd' 0(n - ^/ 3 ) 


(43) 


The second of the series on the right converges uniformly and the series 


sin 2nd sin 2nd’ 1 \ cos 2n(d — i ' ) 1 \ cos 2n(d + d' ) 

^ , n 2/3 ~ 2 n 2/5 2 X_, ^3 

n=l n=-l n=l 

(44) 


may he summed in expl ioit form. 

It is sufficient to consider the series 

CO 

— — / n * 
n-1 

From the known foimua for the F-f unction 

p CO 

/ o - ^ t 2 " 1 dt = T(z) 

d 0 


(15) 


we have 


~iO0 

'lAn- 2 / 3 e niC? = / x- 2 / 3 dx 


ll r 

\ / 


J, 


(46) 


whence 


v~ 


•*\CT 


n il > n- 2 / 3 e nicp = j x-h3 _2'" icp 


n=2 


j 


0 


1 — e 


-x+f 




dx 


(47) 


Formula (47) is obtained from (46) by summing the geometric series 
taking account of the fact that 
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1 _ e " x+iCp | ^ sin q? for 


0 < cp | 

1 < cp < 2jt 


1 


c+iCp l < sin 9 for 


<'f< ? r 


and therefore 



e W( -x+icp) 


1 - 


-x+i<P 

e 


dx <. 


1 

jsin <$> | 




e 


-Wx 


dx = 



jain C P| 



— -> 0 for N — > 0 


or corre spend ingly 


r-l/3 


,N( -x+icp) 
-x+i 1 ? 


dx g: 


JO 



Taking the real part of formula (47) we obtain the required sum of 
aeries (45) : 


fz\ ' cos rfP 

r W zl T?r$ 

n=l 


--1/3 


s ~ x cos cp - e 


-2x 


1 - 2e“ x cos cp -i- © 


JO 


-2x 


; dx (48) 


We proceed to the investigation of the properties of the function (48) 
for cp = 0 and cp = 2n . It is sufficient of course to investigate 
the function (48) near CP = 0. We have: 



28 


NACA TM No, 1155 


- 1/2 \^\ ' cos n r P 

wzl, ~ jj * 

n=l 


■ -1/3 e x coa cp - 1 


} 2x _ 2e x cog cp + 1 


dx + 


* 0 


,-l/3 ( 1 + x) cos CP - 1 

(l + x) 2 - 2(1 + x) COS Cp + 1 


dx + 0(1} = 


r 1 


c -l/3 (l + x) cos cp ~ 1 

(l + x) 2 — 2(1 + x) COS Cp + 1 


dx + O(l) = 


[ x 2/3 dx 5 

I of + cp ^cp 

-'0 JO 


=1/3 


z^dz 
z 6 + 1 


+ 0 ( 1 ) 


(49) 


The last integral is most simply computed with the aid of residues*. 
Further, 


.q>-l/3 


z 4 dz 1 
2 


z 6 + 1 


-'CD' 


rl/3 


z~dz 
z 6 + 1 


l|2«i ^-rti/6 + e 1 /2 + e -5ni/6] + 0 ( C pl/3 ) \ = | + 0 (cpl/3) 

(50) 

(The first terms come from the residues, and the term 0(cp- L / 3 ) from 
the integral over a semicircle of radius r P~V 3 ) 


For this remark which greatly simplifies the preliminary deri- 
vation the author is indebted to A. Nikolsky. 
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Thus ve have finally 


p/2 W ' cos ncp 


%z n - cp 


+ 0 ( 1 ) 


(51) 


whjch gives for the kernal K(0, 0‘ ) 


K(e, 0 1 ) = 


Mo) 


’(|)v(0) 


,1/3 


( 

l 1 


9 - 9 ' I 


- ( 2 9 0 - 9 - 9 ' 


-1/3" 


) 


-1/3 _ 


+ 0 ( 1 ) 


(0 + 9 ' ) 


-1/3 


(52) 


That the obtained kernel K(0 ; 0') actually expresses the 
boundary values of A 1 ( 9 , 0) in terms of ^ a (0 ; 0) is established 
first in the case where A>q is expressed through a finite trigono- 
metric series and in the second case by passing to the limit making 
use of the re spre sentation of in the form of the Chaplygin series. 
Thus the second lemma has been proved. A similar lemma has been 
proven by Tricomi for the equation 


a 2 . 


dx 5y 


s -4-o 


It plays an essential part in reducing the problem considered to an 
integral equation of Fredholm of the second kind. 

We now proceed to the question of the possibility of representing 
the solution of the problem of flow from a vessel with the aid of 
series of the type of Chaplygin. We recall the formulation of this 
problem. There is sought a bounded solution of equation (3) in a 
region of the (u, v) plane (see fig. 4) such that 

= 0 on 0C, \J/ = - | on OA'B' (53) 


We consider now a second solution of equation (3) defined 

by the equation 


Ar 


¥ 2 0 O 


(54) 


This solution satisfies the following conditions of Tricomi: 
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r = 
v = 

if = 


0 on 0C 

0 on 0A* 

Q e __ Q 

2 0 O 2 


> 

on A'B’ 


We assume that this solution exists and that for 
fies the condition proven hy Tricomi in the case 


g(e) 


bo 


0=0 


O(0- 1 / 3 ) 


or in general that the square of the function g 
Then 

CO 

g(e) = / b n sin zve 
n=l 


where 



CO 

and / h 2 n converges. 
n=l 

Then according to Chaplygin the solution 1 
the circle OCA’ will he 




CO 



^n £y(0) 


sin 2U9 


where 

t’p(O) 

An ” Qo) 


In particular on the arc CA 1 


In the case of analytic Laval nozzles this c 
holds. 


(55) 

0 = 0, chiw/cTi' sati 
of equation (])* 

(56) 

9 ) he integrahle . 

(57) 

(57a) 

in the sector of 

(58) 

(53a) 


condition actually 
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^*(o, e) = f(e) = 



2ue 


(59) 


Ab the results obtained by Chaplygin have shown, the convergence of 
the series (58) in the sector OCA* is assured. In the character- 
istic triangle CA'B 1 it is as yet impossible to say anything as 
regards the convergence of the series on the basis of these results. 
We have shown, however, in another paper (reference 7) that for con- 
tinuous variation of the Cauchy data on the arc of the transition 
line (a = 0) the solution of the equation of Tricomi in the corre- 
sponding characteristic triangle varies continuously. From this it 
follows that the problem of Tricomi stated by us is solved in the 
form of the series 


* = ~ ? 



*n 


u°> 


sin 2v6 


( 60 ) 


where the coefficients 

UCo(e)] 


a n are determined from the condition 


CO 

/ a n 
n=l 


>1h 

xyo) 


cm 2ve = | - l) for h 


< 0 < Or 


(61) 


0 = o(0) denoting the dependence of a on 9 along the arc of the 
epicjfdoid . 

Thus out problem under the assumptions made has been reduced to 
the solution of an infinite system of ordinary linear equations. It 
may be attempted to solve this system approximately keeping only a 
finite number of terms in the infinite sum (equation (61)) and requiring 
only a correspondingly finite number of chosen values of d. 


Translation by S. Reiss, 
National Advisory Committee 
for Aeronautics. 
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